Energy conservation in approximated solutions of multiple beam scattering problems 
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Solutions in the form of series expansion, as the Born approximation, are very useful for describing 
time-independent scattering of quantum particles. In this work, it is mathematically demonstred 
that such solutions, when applied to multiple beam scattering phenomena, can lead to energy 
violation where the number of incident and scattered particles is not preserved. General basic 
conditions for developing consistent solutions of time-dependent multiple beam scattering problems 
are outlined. 
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Elastic scattering of quantum particles by extended po- 
tentials is a well-known phenomenon in Physics. It is 
extremely important in Material Science where beams of 
particles, like photons and electrons, are used to probe 
the atomic structure of the matter. Such phenomena are 
describable by scattered waves - vectorial electromag- 
netic waves for photons and scalar wave functions for 
matter particles - written in the form of series expan- 
sions, for instance 



tf(r) = Mr) + Mr) + Mr) 
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as been the probability amplitudes accounting for n th - 
order scattering events. The number of terms is related 
to the extension of the potential, which increases higher- 
order scattering probabilities. Eq. (JIJ stands for the 
general form of the Born series approximation for time- 
independent scattering potentials; and although it is con- 
ceptually simple and of practical usage, it is quite lim- 
ited for describing multiple wave scattering problems. In 
particular, to describe the excitement of one wave when 
others are already excited since such a situation would 
configure time-dependent scattering regime. 

The first objective of this article is to analytically 
demonstrate that phenomenological interpretation of 
multiple wave scattering based on Eq. (JIJ leads to en- 
ergy violation where the potential scatters more parti- 
cles than exist in the incident beam. Consequently, the 
Born series expansion should not be taken as a generally 
valid form of solution for quantitative and accurate de- 
scription of this type of phenomenon. This remark does 
not depend on the number of terms considered in the 
series expansion. The second objective is to develop a 
general time-dependent solution of the phenomenon that 
preserves the number of particles entering and leaving 
the volume occupied by the scattering potential. 

In the standard Born approximation, as find in text 
books 0], ipo(r) stands for the incident wave everywhere. 
It is a reasonable approach for perturbation theories 
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FIG. 1: General configuration of multiple beam elastic scat- 
tering of quantum particles (photons, electrons, neutrons, 
ions,- • • , etc) by extended potentials (macroscopic samples). 
It is characterized by scattered waves, 9g(t), within distinct 
non-overlapping directions of the space. The power of the 
scattered beams, i.e. the total number of particles per unit 
of time, are given by I G = § *G(r)*G( r ) r 2 d9. (G = 0, A, 
B, • ■ ■ ) where the integral is carried out on any closed surface 
outside the potential's volume. 



where the scattered beams are much weaker than the in- 
cident one. However, for a general treatment considering 
very strong scattered beams, ipo(r) has to be the incident 
wave inside the potential's range only; outside, it must 
give the probability amplitude for non-interacting parti- 
cles, those passing through the potential without been 
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scattered. 

A multiple beam scattering configuration, as schema- 
tized in Fig. ^ occurs when Eq. can be rewritten as 

tf(r) = * (r) + *A(r) + * B (r) + --- (2) 

since 

oo 

* G (r) = XI ^n,G(r), G - 0, A, B, • • • , (3) 

and 

^i,fl(r) = ^n(r) (4) 

G 

where ipo,o(r) = V'o( T *)- Since none contribution of the 
incident wave appears in the directions of the scattered 
ones, ipo,G^o( r ) = 0; and neither l st -order scattering in 
the direction of the incident wave, ipifl(r) = 0. Moreover, 
to characterize a multiple beam scattering configuration, 
it is also necessary that 

j ( P )¥*(r) r 2 dr> = l„ + l A + I B + • • • = N (5) 

where 

Iq = j^ G {r)^* G {r) r 2 dfl (6) 

provides the number of particles per unit of time scat- 
tered within distinct non-overlapping solid angles. It 
guarantees that interference of probability amplitudes oc- 
curs exclusively among the series expansion terms of each 
x Pg( t *)i Eq- ©■ The integrals are carried out on any 
closed surface outside the potential's volume. 



The challenge in describing the scattering of multiple 
beams by approximated solutions can be summarized in 
the above example, in how to go from Eq. J2J to Eq. (jHJ 
without violating the equality in Eq. @. In more spe- 
cific words, it is necessary to described not only how the 
extra terms in Eq. (jSJ) are excited, but also how the terms 
already excited in Eq. (7J are affected by the excitement 
of the new beam, Ib in this case. For a qualitative de- 
scription, one may assume that ipQ = ip' , i[>^q = ^2 cm 
and ipi,A = "01 A' an d that the extra terms in Eq. (jHJ) 
are switched on as the potential varies in time. However, 
in this phenomenological description the total number of 



Approximated solution in series expansions are always 
valid, in principle, for time-independent potentials, or 
stationary scattering regime. To use the same form of 
solution for describing the excitement of an extra beam, 
besides those already excited, further developments are 
necessary otherwise the incident and scattered number 
of particles are not preserved in the description of the 
process. This point is demonstrated here by means of 
an example. Let us assume a small potential where 3rd- 
order scattering terms are negligible, as well as higher 
order ones, i.e. VVt>3,G ~ 0. According to Born approxi- 
mation, 

*aW = KaW ( 7 ) 

are the scattered waves under a two-beam excitement 
condition, Iq and l' A . On the other hand, when another 
beam is excited, Ib for instance, the scattered waves are 
given by 

*o(r) = Mr)+^(r) + ^(r) 

* A (r) = Vi.aM+CaW ( 8 ) 

*n(r) = V-lbW + CbM • 
(A,B) superscripts are used on 2 nd -order waves to iden- 
tify the l s *-order ones originating them. For instance, 
■02*0 ( r ) stands for the rescattering of ipi,A(f) towards the 
forward-transmitted wave, ipo( r )- 

The number of particles entering and leaving the ef- 
fective range of the potential must be preserved under 
pure elastic scattering. It means that for an incident 
beam constant in time, No in Eq. JSJ) is constant, and 
therefore, the total scattered intensities of the waves in 
Eqs. and JHJ are the same, i.e. I + Ia = Io + Ia + Ib- 
It leads to 



(9) 

I 

scattered particles is not preserved; unless for a very par- 
ticular coincidence where V^g wou ld provide destructive 
interference effects with the other terms of ^G( r ) (G = 
or A) so that these destructive interference would account 
for the exact number of particles in Ib- Although some 
phase relationships may exist among the scattered waves, 
they vary from one potential to another since the phases 
are intrinsically related to the internal three-dimensional 
structure of the scattering potentials. In the next part of 
this article, a general time-dependent solution, in form of 
series expansion, is proposed for describing the multiple 
beam scattering phenomenon. 



I 

+ Cof + ICaI 2 } r2dn = j {l^o + ^ + Col 2 + hk,A + Cil 2 + IV>i,b + CbI 2 } r2dn ■ 
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Conservation of the number of scattered particles in 
approximated solutions is possible by taken into account 
the scattering probabilities of the particles when traveling 
through the potential. By introducing pn a(n, t) as the 
scattering probability from beam G to beam H after n 
scattering events, the population of particles in the beam 
H that have already been scattered (bounced) from one 
beam to another n + 1 times is 



P H (n+l) 



G 



Ph,g(M)Pg(«-) 



(10) 



where pH.,G(n,t) = for H = G, and the time, t, depen- 
dence of these probabilities are determined by the poten- 
tial, i.e. on how it varies in time until a multiple beam 
scattering configuration is achieved. Moreover, assuming 
a slow time variation to assure that at any time instant 
the populations are given by Eq. I(10|) . the behavior of 
the populations as a function of n can be inferred by the 
sum of probabilities, 



(11) 



When sq(u) = 1 the n-bounced particles in the beam G 
have 100% probability to be scattered towards another 
beam before leaving the potential's range. Therefore, 
outside the potential, the total number of particles on 
beam G is calculated as 



Ig 



Ei 



l-s G (n)]P G (n) 



(12) 



and Sb = Po.b + Pa,b are taken as constant values as a 
function of n. Then, by upgrading the set of scattered 
waves in Eq. (JHJ) according to the new definition given in 
Eq. 1(13(1 . and replacing then into Eq. © we obtain 

I 2t N (l -PA,0 -PBfi +P0,APA,0 +PO,BPB,o) 

I a — N (pa,o + Pa,bPb,o - Po,aPa,o - Pb,aPa,o) 
Ib — N (pb,o + Pb,aPa,o - Po,bPb,o - Pa,bPb,o) 

(14) 

where triple products of the ph,g probabilities are disre- 
garded since ip n >3,G ~ for such a small potential. 

According to the above intensities, Eqs. Q14JI. the total 
number of particles is always preserved, i.e. 



I A + Ib k N 



(15) 



independently of the number of scattered beams. For in- 
stance, when pb,o = Pb,a = beam B is not excited, but 
the sum of intensities still provide the incident number 
of particles, I + I A ~ No. If beam A is also switch off, 
PA,o = Pa,b = 0, no particle-potential interactions occur 
so that I ' = N . 

Strong experimental evidences suggesting scattered 
waves as given in Eq. 1)13(1 are found in recent X-ray mul- 
tiple diffraction experiments Q • There, the linear polar- 
ization of synchrotron radiation is exploited to enhance 
the Aufhellung effect 0, • In terms of scalar wave func- 
tions, such an experiment was equivalent of monitoring 



It implies that when Sg(^) = 1 none n-bounced par- 
ticles are found on beam G. On the other hand, when 
sc(n-) <C 1 the potential scatters a significant fraction of 
the population of n-bounced particles through beam G. 

A direct comparison between Eqs. JJjJl and (|12[1 leads 
to a straightforward conclusion that the ^E , g( t ') waves, as 
given in Eq. J3J), must be redefined outside the potential 
as 



Voir) = ^2 Vl - s gM ipnM r ) 



(13) 



The physical meaning of yl — sc(n) is clear. It re- 
duces the probability amplitudes of scattering towards 
beam G when other beams are excited, i.e. when the par- 
ticles on beam G, inside the potential, have a non null 
probability of leaving the potential's volume via other 
beams. Let us use the three-beam configuration repre- 
sented in Eq. © to demonstrated that the above form of 
solution, Eq. (|T3l . can be used to described the excite- 
ment of the third beam and at the same time preserving 
the total number of scattered particles, No. 

The small potential responsible for the scattered waves 
in Eq. © can be represented by a set of scattering prob- 
abilities, so that s = pa,o + Pb,o, sa = Pq.a + Pb,a, 



(1-sa) 



(B) 

2. A 



r 2 dfl 



(16) 



when exciting Ib- The outcome is a well defined intensity 
reduction owing to the relative enhancement of sa if the 
contribution of V^A ^ s eliminated. 

In summary, multiple wave scattering configurations 
are mathematically defined. It demonstrates that the 
number of incident and scatted particles is not preserved 
when time-independent solutions in form of series expan- 
sion are extended to describe time-dependent scattering 
problems. Basic requirements of time-dependent solu- 
tions are pointed out. It could also be used, in princi- 
ple, to upgrade available solutions to account for the 
Aufhellung effect, commonly observed in multiple beam 
X-ray diffraction. 
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